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Abstract. We present a new theoretical framework for the study of B — > X u lv decays, which 
includes all known perturbative and non-perturbative contributions and a description of leading 
and subleading Fermi motion effects. The perturbative and non-perturbative regimes are separated 
by a "hard" Wilsonian cutoff fi ~ 1 GeV. We bring into focus some problems related to the high q 2 
region and to Weak Annihilation effects. We provide estimates of the CKM parameter \V u i,\ using 
the described framework and discuss the related theoretical uncertainty. 



1. Introduction 

The precise determination of the element \V u b\ of the CKM matrix is an important test of the 
flavour structure of the Standard Model (SM) and is crucial in the indirect search for New Physics. 
The latest global fit to the Unitarity Triangle (UT) including all flavour changing observables but 
a direct determination of \V u b\ predicts \V u b\ = (3.44 ± 0.16) x 10~ 3 [lj. This value agrees within 
errors with the exclusive determination, that relies on lattice QCD or light-cone sum rules [3] 
and that is still affected by somewhat large theoretical errors. An inclusive analysis is in principle 
the cleanest method to precisely determine \V u b\. This is based on the comparison between the 
decay rate of B — > X u lv measured by experiments and the corresponding theoretical prediction. 
The latest HFAG world average 0] yields an inclusive \V u b\ which is about 2.5a higher than the 
value preferred by the global UT fit, calling for a deeper investigation of the process. 

The theoretical description of B — > X u iv is based on a local Operator Product Expansion. 
Inclusive quantities are organized in a double series in a s (perturbative QCD corrections) and in 
l/m& (Heavy Quark Expansion). The very same method was successfully applied to the b — > c 
decay and led to a precise determination of \V c b\, within 2%. The description of charmless decays is 
more involved due to the dominant charmed background that needs to be rejected by experiments 
imposing very stringent cuts. These cuts can spoil the convergence of the OPE introducing 
sensitivity to nonlocal effects, such as the motion of the b quark inside the meson (Fermi motion), 
that can be parameterized in terms of a light-cone distribution function (or "shape function"). 
The lowest integer moments of the distribution function are constrained by the OPE [5] and they 
are expressed in terms of the b quark mass and of the same 5 and 6 dimensional operators that 
contribute to B — > X c lv. Such expressions are universal, i.e. independent of the process, and 
shared by the radiative decay B — > X s ~f only as long as l/m& corrections are neglected. 




An OPE-based treatment of shape function effects in B — > X s ^ including subleading (1/m;,) 
effects was developed in [7] and turned out to describe well experimental data. In [8] a similar 
procedure was undertaken for the case of semileptonic decays, where many complications arise, 
mostly due to the kinematics taking place at different q 2 . In the following we illustrate the main 
features of this procedure and show some meaningful results. 



2. Theoretical framework 

2.1. Perturbative corrections in a Wilsonian approach 

All observables describing the B — > X u £v decay can be extracted via integration over the triple 
differential width: 



where qo and Eg are the total leptonic and the charged lepton energies in the B meson rest frame, 
q 2 is the leptonic invariant mass and W1-3 are the three structure functions relevant in the case of 
massless lepton. 

Perturbative corrections to the structure functions to order 0(a s ) have been known for 

quite long [9], whereas O{a 2 (5o) corrections recently appeared in [ID] . Both calculations were 
performed in the on-shell scheme. 

It has been stressed several times in the literature jTT] that an on-shell definition of the b quark 
mass becomes ambiguous as soon as power suppressed terms are included. The pole mass is better 
traded with a running mass m^). To this purpose, perturbative corrections to order O(a 2 (3o) 
are calculated anew in [8] in the presence of a "hard" Wilsonian cutoff fi ~ 1 GeV. This new 
scale separates the perturbative regime of gluons with energies larger than fi from the "soft" (non- 
perturbative) regime of gluons with energies lower than fj,. The contributions of soft gluons are 
then absorbed into a redefinition of the heavy quark parameters, consistent with the way they are 
extracted from fits to B — > X c tv moments in the kinetic scheme |12| IB]. Physical observables are 
of course independent of the cutoff. 



2.2. Fermi motion 

As already mentioned in the Introduction, Fermi motion is encoded in a distribution function, 
whose lowest integer moments are constrained by the local OPE. As soon as l/mj, corrections 
are retained, such moments cease to be universal: they have different expressions for each of the 
three structure functions in eq. ([I]) and show an explicit q 2 dependence. To preserve generality we 
introduce three separate distribution functions, one for each of the structure functions, depending 
on the light-cone component of the b quark momentum on q 2 and on the Wilsonian cutoff 

(/j,). Hadronic structure functions are then defined via a convolution of the perturbative structure 
functions with the distribution functions, whose expression is derived at the leading order in 1/rrib 
and a s and assumed to be valid also at higher orders: 



Wiiqo, q 2 ) ~ J dk + F l (fc + ,g 2 , / x) [go - ^ (l - ,q 



Model-dependence resides only in the Ansatz employed for the distribution functions. In the analysis 
of |8j a set of about 80 different functional forms, inspired by those already present in the literature 
(exponential, Gaussian, Roman, hyperbolic), is tested and the related uncertainty on \V u b\ turns 
out to be rather small (see Sec. [3|). 
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Table 1. Values of \ V u b \ obtained using different experimental results and their experimental and theoretical 
uncertainties (in percentage) due to various sources (see text). / is the estimated fraction of events. 



2.3. The high q 2 region 

The impact of Fermi motion becomes irrelevant at high q 2 and the developed formalism is no 
more applicable. The OPE itself shows a number of pathological features in this kinematical 
region, due to the emergence of unsuppressed higher order terms. For instance, the OPE predicts 
a value for the variance of the distribution functions which decreases at increasing q 2 , reaching 
even negative values. Moreover, Wilson coefficients of power suppressed operators become more 
and more important and already the coefficient of the Darwin term (p^) shows a divergence at 
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This singularity is removed at the level of the total rate by a one-loop penguin diagram that mixes 
the Weak Annihilation (WA) four-quark operator into the Darwin operator |13|. 114]. However, as 

region 



we are interested in differential distributions as well, a dedicated treatment of the high-g 2 



(q 2 > Q* ~ 11 GeV 2 ) is mandatory. In [8] two different methods are proposed and their difference 
is used to estimate the associated uncertainty: 

a) we model the tail in a way consistent with positivity of the spectra and including a WA 
contribution (X) through a Dirac-<5 localized at the endpoint (default method): 
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b) we extend the Fermi motion description of low q 2 to the high q 2 region, freezing the shape 
function at q 2 = q 2 and using it in the convolution of eq. ([2]) at higher q 2 . 



3. Results and theoretical uncertainties 

We take advantage of some of the latest experimental measurements to extract values of \V u h\ using 
the descripted framework. We leave the task of an average of these results to a future, hopefully 
dedicated, experimental analysis. We consider: 

A Belle analysis with M x < 1.7 GeV and E e > 1.0 GeV [E]; 

B Belle and Babar analyses with M x < 1.7 GeV, q 2 > 8 GeV 2 , and E t > 1.0 GeV [151 US]; 
C Babar analysis with E t > 2.0 GeV [17] 

Results are summarized in Table 1 and were obtained from a C++ code available upon request. 
The reported values of \V u b\ are obtained using the default setting, namely an exponential Ansatz 
for the distribution functions, the prescription a) for the high q 2 tail at X = (see previous section) 
and the central values of the fit in [6] as input parameters, at fx = 1 GeV. / is the estimated fraction 
of events, whereas the following columns show different sources of uncertainty, namely: 



• Experimental error (exp). 

• Parameteric error (par): it is extracted taking into account all correlations between non- 
perturbative parameters [6] and varying a s = 0.22 by ±0.02 as uncorrelated. The uncertainty 
on mf, is by far dominating. 

• Perturbative error (pert). 

• Errors related to the high q 2 region: we consider the difference between methods a) and b) 
(tail model), the value of q 2 at which the modelling sets in (q 2 ) and the error due to WA 
effects (X). We let X vary in a range consistent with the 90% confidence level bound set by 
CLEO on the size of WA [15] : namely < X < 0.04. We stress that the error related to X is 
asymmetric and points to a lower value of \V u b\. 

• Functional form dependence (ff), estimated using about 80 different Ansatze for the distribution 
functions. 

It is worth stressing that all combinations of cuts considered include the high q 2 region discussed in 
the previous section which, as we showed, is plagued by poorly controlled effects. However, Belle 
measurements A and B can be easily combined in order to obtain an estimate of \ V u b\ with an upper 
cut on q 2 , namely for the combination Mx < 1-7 GeV, Et > 1.0 GeV, and q 2 < 8 GeV 2 . This yields, 
within our framework, a value of \V u b\ much lower than in all other cases (\V u b\ = 3.18 x 10 -3 ). 
and it might signal some bias in the treatment of the high q 2 region either on the experimental 
or on the theoretical side. Only a dedicated experimental analysis with an upper cut on q 2 could 
probably shed some light on this issue. 



4. Summary and References 

We presented a new approach for dealing with the triple differential width of B — > X u lu decays, in 
a framework characterised by a hard Wilsonian cutoff [i ~ 1 GeV. The method developed takes into 
account all known perturbative and non-perturbative corrections. Fermi motion is treated at the 
subleading level as well. Some problems related to the high q 2 region of the process were pointed 
out, that were probably underestimated in the past and that still deserve a deeper investigation. 
We also presented some numerical results with the associated uncertainties and put forward the 
suggestion of a new experimental analysis with an upper cut on q 2 . 
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